Single-Molecule Investigation of the Binding Interface Stability of SARS-CoV-2 Variants with ACE2.

The SARS-CoV-2 pandemic spurred numerous research endeavors to comprehend the virus and mitigate its global severity. Understanding the binding interface between the virus and human receptors is pivotal to these efforts and paramount to curbing infection and transmission. Here we employ atomic force microscopy and steered molecular dynamics simulation to explore SARS-CoV-2 receptor binding domain (RBD) variants and angiotensin-converting enzyme 2 (ACE2), examining the impact of mutations at key residues upon binding affinity. Our results show that the Omicron and Delta variants possess strengthened binding affinity in comparison to the Mu variant. Further, using sera from individuals either vaccinated or with acquired immunity following Delta strain infection, we assess the impact of immunity upon variant RBD/ACE2 complex formation. Single-molecule force spectroscopy analysis suggests that vaccination before infection may provide stronger protection across variants. These results underscore the need to monitor antigenic changes in order to continue developing innovative and effective SARS-CoV-2 abrogation strategies.

Molecular Insight into the Self-Assembly Process of Cellulose Iß Microfibril.

The self-assembly process of ß-D-glucose oligomers on the surface of cellulose Iß microfibril involves crystallization, and this process is analyzed herein, in terms of the length and flexibility of the oligomer chain, by means of molecular dynamics (MD) simulations. The characterization of this process involves the structural relaxation of the oligomer, the recognition of the cellulose I microfibril, and the formation of several hydrogen bonds (HBs). This process is monitored on the basis of the changes in non-bonded energies and the interaction with hydrophilic and hydrophobic crystal faces. The oligomer length is considered a parameter for capturing insight into the energy landscape and its stability in the bound form with the cellulose I microfibril. We notice that the oligomer-microfibril complexes are more stable by increasing the number of hydrogen bond interactions, which is consistent with a gain in electrostatic energy. Our studies highlight the interaction with hydrophilic crystal planes on the microfibril and the acceptor role of the flexible oligomers in HB formation. In addition, we study by MD simulation the interaction between a protofibril and the cellulose I microfibril in solution. In this case, the main interaction consists of the formation of hydrogen bonds between hydrophilic faces, and those HBs involve donor groups in the protofibril.

Assessing the Stability of Biological Fibrils by Molecular-Scale Simulations.

The nanomechanical characterization of several biological fibrils that are the result of protein aggregation via molecular dynamics simulation is nowadays feasible, and together with atomic force microscopy experiments has widened our understanding of the forces in the regime of pN-nN and system sizes of about hundreds of nanometers. Several methodologies have been developed to achieve this target, and they range from the atomistic representation via molecular force fields to coarse-grained strategies that provide comparable results with experiments in a systematic way. In this chapter, we discuss several methodologies for the calculation of mechanical parameters, such as the elastic constants of relevant biological systems. They are presented together with details about parameterization and current limitations. Then, we discuss some of the applications of such methodologies for the description of bacterial filament and ß-amyloid systems. Finally, the latest lines of development are discussed.

Martini 3 Model of Cellulose Microfibrils: On the Route to Capture Large Conformational Changes of Polysaccharides.

High resolution data from all-atom molecular simulations is used to parameterize a Martini 3 coarse-grained (CG) model of cellulose I allomorphs and cellulose type-II fibrils. In this case, elementary molecules are represented by four effective beads centred in the positions of O2, O3, C6, and O6 atoms in the D-glucose cellulose subunit. Non-bonded interactions between CG beads are tuned according to a low statistical criterion of structural deviation using the Martini 3 type of interactions and are capable of being indistinguishable for all studied cases. To maintain the crystalline structure of each single cellulose chain in the microfibrils, elastic potentials are employed to retain the ribbon-like structure in each chain. We find that our model is capable of describing different fibril-twist angles associated with each type of cellulose fibril in close agreement with atomistic simulation. Furthermore, our CG model poses a very small deviation from the native-like structure, making it appropriate to capture large conformational changes such as those that occur during the self-assembly process. We expect to provide a computational model suitable for several new applications such as cellulose self-assembly in different aqueous solutions and the thermal treatment of fibrils of great importance in bioindustrial applications.

Mapping Mechanostable Pulling Geometries of a Therapeutic Anticalin/CTLA-4 Protein Complex.

We used single-molecule AFM force spectroscopy (AFM-SMFS) in combination with click chemistry to mechanically dissociate anticalin, a non-antibody protein binding scaffold, from its target (CTLA-4), by pulling from eight different anchor residues. We found that pulling on the anticalin from residue 60 or 87 resulted in significantly higher rupture forces and a decrease in koff by 2-3 orders of magnitude over a force range of 50-200 pN. Five of the six internal anchor points gave rise to complexes significantly more stable than N- or C-terminal anchor points, rupturing at up to 250 pN at loading rates of 0.1-10 nN s-1. Anisotropic network modeling and molecular dynamics simulations helped to explain the geometric dependency of mechanostability. These results demonstrate that optimization of attachment residue position on therapeutic binding scaffolds can provide large improvements in binding strength, allowing for mechanical affinity maturation under shear stress without mutation of binding interface residues.

Molecular insights into receptor binding energetics and neutralization of SARS-CoV-2 variants.

Despite an unprecedented global gain in knowledge since the emergence of SARS-CoV-2, almost all mechanistic knowledge related to the molecular and cellular details of viral replication, pathology and virulence has been generated using early prototypic isolates of SARS-CoV-2. Here, using atomic force microscopy and molecular dynamics, we investigated how these mutations quantitatively affected the kinetic, thermodynamic and structural properties of RBD-ACE2 complex formation. We observed for several variants of concern a significant increase in the RBD-ACE2 complex stability. While the N501Y and E484Q mutations are particularly important for the greater stability, the N501Y mutation is unlikely to significantly affect antibody neutralization. This work provides unprecedented atomistic detail on the binding of SARS-CoV-2 variants and provides insight into the impact of viral mutations on infection-induced immunity.

Characterization of Structural and Energetic Differences between Conformations of the SARS-CoV-2 Spike Protein.

The novel coronavirus disease 2019 (COVID-19) pandemic has disrupted modern societies and their economies. The resurgence in COVID-19 cases as part of the second wave is observed across Europe and the Americas. The scientific response has enabled a complete structural characterization of the Severe Acute Respiratory Syndrome-novel Coronavirus 2 (SARS-CoV-2). Among the most relevant proteins required by the novel coronavirus to facilitate the cell entry mechanism is the spike protein. This protein possesses a receptor-binding domain (RBD) that binds the cellular angiotensin-converting enzyme 2 (ACE2) and then triggers the fusion of viral and host cell membranes. In this regard, a comprehensive characterization of the structural stability of the spike protein is a crucial step to find new therapeutics to interrupt the process of recognition. On the other hand, it has been suggested that the participation of more than one RBD is a possible mechanism to enhance cell entry. Here, we discuss the protein structural stability based on the computational determination of the dynamic contact map and the energetic difference of the spike protein conformations via the mapping of the hydration free energy by the Poisson-Boltzmann method. We expect our result to foster the discussion of the number of RBD involved during recognition and the repurposing of new drugs to disable the recognition by discovering new hotspots for drug targets apart from the flexible loop in the RBD that binds the ACE2.

Quantitative determination of mechanical stability in the novel coronavirus spike protein.

We report on the novel observation about the gain in nanomechanical stability of the SARS-CoV-2 (CoV2) spike (S) protein in comparison with SARS-CoV from 2002 (CoV1). Our findings have several biological implications in the subfamily of coronaviruses, as they suggest that the receptor binding domain (RBD) (∼200 amino acids) plays a fundamental role as a damping element of the massive viral particle's motion prior to cell-recognition, while also facilitating viral attachment, fusion and entry. The mechanical stability via pulling of the RBD is 250 pN and 200 pN for CoV2 and CoV1 respectively, and the additional stability observed for CoV2 (∼50 pN) might play a role in the increasing spread of COVID-19.

Energy Gaps of Polyradicals from an Effective and Transferable Hamiltonian with through-Bond Interactions.

Current model Hamiltonians and ab initio many-body quantum treatments of π-conjugated polyradicals formed from hydrocarbons produce divergent results because of numerical complexity and large size of the basis-function set used. We propose an alternative, three-term Hamiltonian, to describe these various polyradicals that simplifies considerably the computational cost while providing a physical interpretation for all three terms and a high degree of model universality. The essential feature of this Hamiltonian is a term, not present in previous models, describing the three-sited through-bond interaction that governs the noninteracting spin-up and spin-down sectors. A computation of the lowest energy gaps and spin configurations for the smaller polyradicals demonstrates the efficacy of the model and its potential in applications in revealing electrical conductivity and ferromagnetism of the more complicated substituted polyradicals.

On the separability of the extended molecule: Constructing the best localized molecular orbitals for an organic molecule bridging two model electrodes.

Based on a quantum chemical valence formalism that allows the rigorous construction of best-localized molecular orbitals on specific parts of an extended system, we examined the separability of individual components of model systems relevant to the description of electron transport in molecular devices. We started by examining how to construct the maximally localized electronic density at the tip of a realistic model of a gold electrode. By varying the number of gold atoms included in the local region where to project the total electronic density, we quantitatively assess how many molecular orbitals are entirely localized in that region. We then considered a 1,4-benzene-di-thiol molecule connected to two model gold electrodes and examined how to localize the electronic density of the total system in the extended molecule, a fractional entity comprising the organic molecule plus an increasing number of the closest metal atoms. We were able to identify in a rigorous manner the existence of three physically different electronic populations, each one corresponding to a distinct set of molecular orbitals. First, there are those entirely localized in the extended molecule, then there is a second group of those completely distributed in the gold atoms external to that region, and, finally, there are those delocalized over the entire system. This latter group can be associated to the shared electronic population between the extended molecule and the rest of the system. We suggest that the treatment here presented could be useful in the theoretical analysis of the electronic transport in nanodevices whenever the use of localized molecular states are required by the physics of the specific problem, such as in cases of weak coupling and super-exchange limits.